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ABSTRACT 


A nonlinear maneuvering model based on ship geometric and mass properties is 
developed. The model can be utilized to evaluate maneuvering performance early in 
the design phase. This model is also used in this thesis as a benchmark in order to 
evaluate the accuracy of the simpler Nomoto's model. The latter is faster to simulate 
and is ideally suited for visual simulation studies. Results comparing the relative 
accuracy and speed of implementation of the two models are presented for different 
inputs and geometric properties. An improvement to Nomoto's model is suggested 
which greatly increases accuracy while maintaining high speed of real time 
implementation. Furthermore, a series of parametric studies is performed in order to 
_ evaluate the sensitivity of fundamental maneuvering properties in terms of basic ship 


geometric quantities. 
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I. INTRODUCTION 


Today's military faces many challenges to do more with less. Since the end of the 
Cold War, the topic of discussion on Capitol Hill has been the shrinking military budget. 
However, the military's requirement to maintain full combat readiness has not changed. In 
order to maintain training at the required level, alternatives to costly full scale maneuvers 
must be developed. | 

Current war gaming models use mostly static or simplified dynamic descriptions of 
both attacking weapons and targets. Clearly, this provides an unrealistic description of the 
situation, with respect to the vulnerability of the target to attack, and does not allow for the 
incorporation of standard tactical maneuvers of the target to avoid incoming weapons. What 
needs to be done is to increase the accuracy of current models by incorporation of targets in 
a dynamic, rather than a static manner. Developing an accurate simulation model, will lead 
to ship designers being able to examine the vulnerability of a ship, in real time through the 
simulation model. 

In this thesis, a nonlinear model based on ship geometric and mass properties will be 
developed. The model can be utilized to evaluate maneuvering performance early in the 
design phase. This model will also serve as a benchmark in order to evaluate the relative 
accuracy of a simpler linear model. Chapter II will show the development of the 
mathematical models. Chapter III will discuss the assessment of the models and compare 
them for visual simulation suitability. Chapter IV will evaluate the sensitivity of 
fundamental maneuvering properties in terms of basic ship geometric quantities. Chapter 
V will assess two variations to the linear model, and compare them to the nonlinear model. 


Finally, Chapter VI will discuss conclusions and recommendations for further work. 











Il. MATHEMATICAL MODELING 


In this dynamical model, the planar motion of a ship on an undisturbed free-surface 
of an incompressible fluid will be considered. The basic dynamics of maneuvering will be 
described and analyzed using Newton's equations of motion. 

We will define two reference frames, one set of axes fixed relative to the Earth 
aligned with directions North, East, and Down (denoted by capital letters), and a second set 
of axes fixed relative to the ship, aligned with the longitudinal, the starboard direction, and 


down (denoted by lower case letters). 











Figure 1 Designation of reference frames 


The heading (yaw) angle, refers to the direction of the ship's longitudinal axis, with 
respect to one of the fixed axes. The drift angle refers to the difference between the heading 
and the actual course. Let the instantaneous speed of the ship be denoted by U, then the 


surge and sway velocities are given by: 


u-UcosB , -v--UsinB 





The dynamical model used in this work has been proposed [Ref. 1] and is given by: 


X'-(m'.m'(L/UUcosB/U-BsinB)+(m ‘sm, Yr'sinB 
Y'--(mm'sm(L/UUsinB/U-B cosB )(m ‘n)r ‘cosB 
N'= (+i,,(LIU 2\(Ur UF 


where, 
B is the drift angle 
r is the rotational or yaw velocity 
L,,and 1, are the mass moments of inertia in both ref frames 
m, m,, and m, are the ship's mass, added mass of the x-axis, and the 


added mass of the y-axis, respectively 


All terms of these equations, and the equations of motion in the remaining report, have been 
nondimensionalized using the ship's length L, the ship's draft d, the fluids density p, and the 
ship's speed U, as shown in the Appendix The terms X, Y, and N are the extemal force along 
the x-axis, the external force along the y-axis, and the yaw moment about the ship's center of 
gravity, respectively. 
It is assumed that X, Y, and N have the following decomposition, 
X=X tk prX p 
Y-¥ ¥ pt] p 


N-N ypN pN p 


where the subscripts represent respective contribution due to, 


H - ship's hull 





P - ship's propeller 

R - ship's rudder 
We can assume that during most ordinary maneuvers, the forward velocity is kept constant 
by the propulsion control. Also, in looking at the ship's response to rudder deflection, it is 
assumed that no significant surge accelerations will take place. From these assumptions we 
are left with, 


Y-Y,+Y, 


N-N,N, 


The components of the rudder forces are taken to be of the form, 


Y=--(1+a,)F,cos6 


N,=-(%p+4,X,)F ycos6 


F,--(A,/L@)C,U,2sin(8-(B-2x,7)) 


C,-6.13K,/(Kg*2.25) 


where, 


F,- force acting on rudder 

Ap- rudder area 

K,- rudder aspect ratio 

U,- effective rudder inflow speed (taken to be equal to U) 

Xp- distance between COG of ship and center of lateral force 

X,,- distance between COG of ship and center of additional lateral 
force 


a,,- ratio of additional lateral force 





Ap Kp Xp, Xy, and a, are constants obtained through model tests. 


The lateral force and yaw moment acting on the ship's hull are expressed as follows: 


¥,-¥,B+¥7+¥,,B1B|+¥_7ir|-(Y,,B+¥,,) Br 


N-N, BN rN B | B [Nr |r| +N, 5, BNF) Br 


The lateral force and yaw moment hydrodynamic coefficients are given as, 


Y -m,-1.5sC, 


1 
a rae 


k 
Yyq-2.97-(1-Cy) 


k | 
Ygq*1.25-—(1-C,)+0.5 


¥,-0.17£C,-0.07 
S 


k 
¥_5,-0-75-C,-0.65 


N--0.54k+k? 


Nik 
N..=-0.48~(1-C.):0.66 
4p" 0.48 (1-C,).0. 


N_-0.5sC,-0.09 


N,,--0.25~C,+0.05 
AY 


Nyo,7-57.5(SC,)*+18.4sC,-1.6 





where, 
s-B/L, Slenderness ratio 


k-2d/L, Aspect ratio 


B 


C ee Block coefficient 
LdB 


m,-2sC,,(1+k), Mass coefficient 


The above model may be simplified by linearizing the sway and yaw equations, as done in 


[Ref. 3], to obtain, 


(m-Y )v{(Y,-mx )P=Y v+(¥ -mU)r+Y,6 


(I,-N,)P-(N,-mx,)¥-N v(N-mx UPN, 8 


where 6 is the rudder deflection angle measured according to the convention of Figure (1); 
positive rudder corresponds to a turn to port. These equations of motion can be manipulated 


to produce a pair of decoupled second order equations in r and v, 


TT F(T ,+T,)Fer-K8+KT,6 
7,T,0+(T,+T, ¥+v-K,8+K,T,6 


We will assume that there 1s no appreciable side slip, and hence, the second equation may 
be ignored. What remains 1s known as Nomoto's equation. The two indicies, K and T, are 
related to the hydrodynamic coefficients by the following, 


_(-m)N,-mx UN LY, | 
'  ¥(N-mx, U), (¥,-mU) 





mx NN, mx, ¥-mU) 


*YN,-mxU)-N(F,-mU) 


T - (NV. y MX) P rae 6 ym), 
NY,-YN, 


NY,-YN, 
YL Wmmx,U)-N(¥,-mU) 


where, Y -N-0.0 


Y =m,-1.5sC, 


Y--0.5nk-1.4sC, 


N-L-I 


N a7 (Xt XP N 


N.--0.54k-k? 


1 


m=2sC,(1+4,), i=1,2 9 A,=k (1-5) 





I(2)sC,(1-2 AE C1-1.6s) 
5 


Nomoto's equation expresses the relationship between the ship's turning rate and the rudder 


angle through a second order transfer function, 


r K+KT,s 
6 TT7,s?(T,+T,)s+1 


which can be simplified to a first order transfer function, 





This equation in differential form is known as Nomoto's First Order Equation, 


Trs+r-K6 


where, 


vy be ey 


By using basic principles of ordinary differential equations, the solution for r, in terms of K 


and T, is given by, 


r - K8&(1-e*") 


Having developed these two models, nonlinear and liner, the two will now be compared to 





see which will give a faster calculation time, and which is more accurate for graphical 


Simulation. 
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Ill. MODEL ASSESSMENT 


Having developed the nonlinear and linear (Nomoto's) models, the two were 
compared for speed of implementation and the relative accuracy for different inputs and 
geometric properties. 

The principal dimensions of length, beam, and draft were varied independently and 
the models were also compared utilizing different rudder angle inputs. 

In every instance of the comparison between these two models, the linear model 
proved to be approximately six times faster than the nonlinear model. This fact makes 
Nomoto's model ideally suited for visual simulation studies. Visual simulations may be 
carried out with minimal threat of reduced frame speed, due to calculation processing time. 

A graph of the advance vs. transfer (i.e., a geographical x-y plot of the ship's path) 
is shown in Figure 2 for both linear and nonlinear models and for different rudder angles. 
Advance and transfer are shown in dimensionless form, with respect to the ship length. 
Rudder angle, indicated by d on the graphs, is in degrees. Solid curves correspond to the 
nonlinear simulation, while dotted curves correspond to Nomoto's first order model. It can 
be seen that the latter accurately depicts the ship's path for relatively large rudder angles. 
However, as the rudder angle is decreased the difference between linear and nonlinear 
simulations is likewise more pronounced. This deviation may be attributed to the increased 
side-slip velocities that are encountered during smaller rudder deflections and are neglected 
in Nomoto's model. This is explored in the next chapter. | | 

The results for the yaw rate r are shown in Figure 3. Both yaw rate and time are 
shown in dimensionless quantities, with respect to ship speed and length. It seems that 
higher rudder angles result in slight overshoot in the yaw rate for the nonlinear model. This 
means that for high rudder angles, the yaw rate 1s basically described by a second order 


system, whose damping ratio becomes smaller. 


1] 





The sensitivity of these results for small variations in the ship's draft is demonstrated by 
Figures 4 through 7. A plot of the advance vs. transfer for different values of the draft based 
on the nonlinear model is shown in Figure 4. It appears that small drafts correspond to more 
maneuverable ships with tighter turning circles. The same result can be drawn from the 
linear simulation results of Figure 5. From the yaw rate simulations of Figures 6 and 7, we 
can see that smaller drafts correspond to higher turning rates at steady state, as well as, at 
smaller time constants. The nonlinear simulations also suggest that a smaller draft results 
in smaller damping ratios, as evidenced by the overshoot in the yaw rate time history. 

The sensitivity, with regards to ship beam, is demonstrated by Figures 8 throug: 11. 
In general, it appears that the beam has very little effect on maneuverability. The nonlinear 
model suggests that maneuverability is increasing for increasing beam, although the linear 
model predicts the opposite. It also appears that the yaw rate damping ratio is decreased for 
decreasing values of the ship's beam. 

The sensitivity, with regards to the ship length, is demonstrated by Figures 12 through 
18. In general, both linear and nonlinear models suggest that maneuverability is increased 
with increasing length. This results in tighter turning circles, both in dimensional and 


dimensionless quantities. 
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Figure 2 Rudder Angle comparison between Nonlinear Model and 
Nomoto's Model 
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Figure 3 Yaw Rate comparison between Nonlinear Model and Nomoto's 
Model 
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Transfer 


Figure ¢ Advance vs. Transfer for Nonlinear Model at different drafts 
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Figure 5 Advance vs. Transfer for Nomoto's Model at different drafts 
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Figure 7 Yaw Rate vs. Time for Nomoto's Model at varying drafts 
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Figure 8 Advance vs. Transfer for Nonlinear Model at varying beams 
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Figure 9 Advance vs. Transfer for Nomoto's Model at varying beams 
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Figure 11 Yaw Rate vs. Time for Nomoto's Model at varying beams 
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Figure 13 Dimensional Advance vs. Transfer for Nonlinear Model at 
varying lengths 
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Figure 15 Yaw Rate vs. Time for Nomoto's Model at varying lengths 
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Figure 17 Dimensional Yaw Rate vs. Time for Nonlinear Model at 
varying lengths 
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IV. LINEAR MODEL VARIATION ASSESSMENT 


In Chapter I, the assumption of negligible side slip velocity was made to obtain 
Nomoto's equation. If this assumption was not made, we would have the pair of decoupled 


second order equations from the linearized sway and yaw equations, 


T.T,0+(T,+T,v+v-K 6+K 7,5 


TT, F(T ,+T,)r+r-K8+KT,6 


The relationship between the indicies, in the yaw rate second order equation and the 
hydrodynamic coefficients was shown in Chapter II. The indices of the side slip velocity 


second order equation are related to the hydrodynamic coefficients by the following, 


_ CNN, -mx WY, 
"YW -mx UN, (Y,-mU) 


. (N,- Po) ¢ a (¥,-mx,U)N, 
‘ (N,,- GUY y (Y,-mU)N, 


The second decoupled second order equation led to formulation of Nomoto's equation, and 
the first expresses the relationship between the side slip velocity and the rudder angle through 
a second order transfer function, 

vy KK Ts 
6 T_T,s*(T,+T,)s+1 


which can be simplified to a first order transfer function, 





we 
6 T s+} 
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This equation in differential form reduces to, 
T v+v-K 6 


where, 


GPs EG Be 


Once again, by using basic principles of ordinary differential equations the solution for the 


side slip velocity is, 


t 


— 


v-K &(1-e ” ) 


In Figure 18, a comparison is shown between Nomoto's model and a linear model 
containing the affects of side slip. The two graphs are almost identical, which leads to the 
conclusion that the original assumption of negligible side slip velocity was correct. 

Now we will compare the nonlinear model with the side slip model. Figures 19 and 
20 lead to the same conclusions as Figures 2 and 3, that relatively speaking, the nonlinear 
model is more accurate than the linear model. Figure 21 shows the comparison between the 
Side slip velocity obtained in the nonlinear vs. the linear case. Once again, it is clear that the 
nonlinear model simulates the situation more accurately. 

The side slip assessment leads to the conclusion that the difference between the linear 
and nonlinear models lies not with the side slip velocity, but with the yaw rate. It can be seen 
in Figure 22 that the difference in the error between the linear and nonlinear models is a 
function of the rudder angle. 

By plotting the rudder angle vs. the steady-state yaw rate ratio (nonlinear over linear), 


we obtain a straight line as shown in Figure 23. By calculating the slope and intercept of this 


ZZ 





line and applying it as a correction factor, we obtain a new relationship for Nomoto's 
equation, 


Tr + r= (a-B|6))K6 


By applying this "correction" factor, it can be seen in Figure 24 that we obtain a mode] that 
is more accurate than both Nomoto's model and the side slip model. On the average the 
corrected model is approximately 10% slower than Nomoto's model. 

In Figure 25, it can be seen that the relative error in the yaw rate between the 
nonlinear model and Nomoto's model has been removed in the corrected model. By 
comparison, the conclusion that the corrected model is the most accurate of the models 
presented can be made. Due to the relative accuracy of the corrected model, it is better 
suited to simulate the turning rate while maintaining an acceptable frame speed during visual 


simulations. 
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Figure 18 Advance vs. Transfer comparison between Nomoto's Model and 
Side Slip Model 
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Figure 19 Advance vs. Transfer comparison between Nonlinear Model and 
Side Slip Model 
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Figure 20 Yaw Rate vs. Time comparison between Nonlinear Model and 
Side Slip Model 
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Figure 21 Side Slip Velocity comparison between Nonlinear Model and 
Side Slip Model 
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Figure 22 Yaw Rate vs. Time comparison between Nonlinear Model and 
Nomoto's Model at varying rudder angles 
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Figure 23 Yaw Rate Ratio vs. Rudder Angle plot used to obtain correction 
factor for Nomoto's Model 
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Figure 24 Advance vs. Transfer plot comparing Nonlinear, Nomoto, and 
corrected Nomoto Models 
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Figure 25 Yaw Rate vs. Time plot comparing Nonlinear, Nomoto, and 
corrected Nomoto Models 
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V. PARAMETRIC STUDY 


During the preliminary ship design phase, it is often possible to introduce small 
changes in fundamental hull properties so that maneuverability is improved. Of course, such 
changes may only be instituted after taking due consideration of all aspects and requirements 
of the design problem. | 

In this Chapter, the maneuvering properties of a ship are explored by utilizing the 
nonlinear model described in the previous chapters, since it reflects reality better than the 
linear model. In order to establish a systematic series of runs, the ship's displacement was 
maintained constant. The dispalcement is usually determined by payload considerations and 
it is unlikely that a designer would have the freedom to change it based on maneuvering 
requirements alone. 

Figures 26 and 27 show results that are obtained by keeping both the displacement 
and the beam constant. The length and draft are then changed within plus/minus 5% of their 
nominal values. It can be seen that the longer and shallower hull form has superior 
maneuvering characteristics compared to shorter and deeper draft ships. 

Figures 28 and 29 show results that are obtained by keeping both the displacement 
and the length constant. The beam and draft are then changed within plus/minus 5% of their 
nominal values.. It can be seen that wide, shallow hull forms have superior maneuvering 
characteristics compared to narrow and deeper draft ships. 

Figures 30 and 31 show results that are obtained by keeping both the displacement 
and the draft constant. The length and beam are then changed within plus/minus 5% of their 
nominal values. It can be seen that the long and narrow hull forms have superior 
maneuvering characteristics compared to shorter and wider draft ships. 


Maneuvering performance is increased by: 


- a longer length 
- jess draft 
- a narrower beam for constant draft 


- a wider beam for constant length 
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Therefore, beam does not have a consistent trend and since it is usually set by static 
stability considerations, it should not be altered based on maneuvering alone. The length is 
the most costly dimension to increase. Not only are material cost increased, but dry dock 
size, and shipbuilding considerations must also come into play. The draft on the othe hand 


will have it's own set of considerations such as harbor and canal restrictions. 
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Figure 27 Yaw Rate vs. Time comparison with constant beam 
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Figure 28 Advance vs. Transfer comparison with constant length 
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Figure 29 Yaw Rate vs. Time comparison with constant length 


32 





Transfer 


Figure 30 Advance vs. Transfer comparison with constant draft 
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Figure 31 Yaw Rate vs. Time comparison with constant draft 
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VI. CONCLUSIONS 


In this thesis two maneuvering models were developed and compared. A nonlinear 
model based on ship geometric and mass properties was developed from [Ref. 1]. The 
second, linear, model was developed by linearizing the sway and yaw equations of the 
nonlinear model, as done in [Ref. 3]. The models were compared for speed of 
implementation and the relative accuracy for different inputs and geometric properties. It 
was concluded that though the linear model was approximately six times faster than the 
nonlinear model, it was not accurate enough. 

At this point a variation was made to the linear model by adding the effects of the 
side slip velocity to possibly explain the difference in the accuracy. Upon further inspection, 
it was determined that the difference in accuracy did not lie with the side slip velocity, but 
with the difference in the yaw rates. | 

By plotting the rudder angle vs. the steady-state yaw rate ratio, a correction factor was 
obtained and applied to the linear model ;causing the steady-state error between the yaw rates 
of the linear and nonlinear models to be removed altogether. This "corrected" linear model 
also proved to be relatively more accurate than either the original linear model and the 
nonlinear model. Though this model is 10% slower than the original linear model, it is still 
much faster than the nonlinear model. 

A series of parametric studies was also performed on the nonlinear model in order 
to evaluate the sensitivity of fundamental maneuvering properties in terms of basic ship 
geometric quantities. It was concluded that maneuvering performance is increased by a 
longer length or by less draft. The beam did not show a consistent trend, and is usually set 
by static stability considerations and probably not alterable for maneuvering alone. 

For further work, it is recommended that these models be incorporated into the 
Computer Science Department battlespace simulation system, NPSNET. Initially this system 
was geared towards ground forces and land based conflicts, but recently a naval component 
has been added in an attempt of attaining a more realistic joint services simulator. The 
models should be checked for accuracy, realism, and speed of implementation. Once a 


baseline is obtained, improvements in the best model for visual simulation can be made. 
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APPENDIX. MATLAB CODE 


Nondimensional quantities, denoted by primes, are defined as follows: 


m'=m/(0.5pL7d) 
m,'=m,/(0.5pL7d) 
m,'=m,/(0.5pL*d) 

7 =L,/(0.5pL“d) 
i, =i, /(0.5pL*d) 

'=X/(0.5pLdU?) 

Y'=¥/(0.5pLdU?) 
N'=N/(0.5pL2dU?) 
r=rL/U 
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% COMPARISON OF LINEAR AND NONLINEAR MODELS 
% ON SAME GRAPH(including Nomoto correction & side slip model) 


% 

$ Define parameters of the simulation (Metric units) 

% 

% 

time = 90.0; % Final time of simulation 
at = .1¢ * Sample time of difference eqns 
2 = 3257400. + % Moment of inertia 

Los 250.03 %$ ship’s length 

RHO: =.1025; % Density of salt water 

XG = 10.0; % Center of gravity 

Mass = 147241915.0; $ Ship mass 

m = Mass/(0.5*RHO*L%3); $ Dimensionless ship mass 
izz = Iz/ (0.5*RHO*L*S); %$ Dimensionless moment of inertia 
Xg = XG/L; %* Dimensionless center of gravity 
d= -0.2618; % rudder angle 

B = 40.77; %* ship's breath(twice the beam) 
D = 20.00; $ ship’s draft 

s = B/L; % slenderness ratio 

V = Mass/RHO; % ship’s displacement 

Cb = V/(L*D*B); * block coeff 

k= 2" D/i: %$ hull aspect ratio 

Kr ae % rudders aspect ratio 

AEDS 7420 %$ rudder area 

Di: Gyo % propeller diameter 

eo Spb 

iV.s> 3000¢ % propeller revolution (rps) 
Ctp: = 1.0: 

aly = 0.633*Cb - 0.153; 

ims = -156.2* (Cb*B/L)*2 + 41.6* (Cb*B/L) ~ 1.76; 

wy 2S: ae Ch- ee DS) De 

“0 = 0.6*wp; %* thrust deduction coefficient 
we = 1] - (eps*(l-wp)); 

K = 0.6* (l-wp)/(l-wr); 

tr .= 0.28*Cb -+° 0.255; 

eta = Dp/(sqrt(AR*Kr) ); 

Cn = (6.135 "hry (Kr 2 25) 

lambda6é = (k/s)*(1-1.6%*s); 

I = (1/8) *s*Cb* (1+lambda6) ; 

lambda2 = (k/s)*(1-s/2); 

m2 = 2*s*Cb* (1+lambdaZ2) ; 

lambdal = k; 

ml = 2*s*Cb* (1+lambdal); 

PI = 3.14: 

XR = -Q0.5; 

XH = =0..52 

mi = 2*s*Cb* (1+lambdal) ; 

gamma = -22.2* (Cb*B/L)*2 + 0.02* (Cb*B/L) + 0.68; 

Ue 1. Os 

UR = U; 

% 

%$ LINEAR MODELS 

% 

beta = .25; % drift angle 

Fn= -(AR/(L*D)) *Cn*UR%2*cos(- gamma*beta); % Force acting on rudder 
% Define hydrodynamic coefficients 

% 

Yvl = -0.5*PI*k + 1.4*sS*Cb; 

Nvi = -k; 

Yril. = ml - 1.5*s*Cb; 

Yd1l = -(1+aH)*Fn; 

Nrl = ~-0.54*k + k‘*2; 

Ndi = -(xR + aH*xH) *Fn; 

Nridot = Izz-I; 

Yvidot = m-m2; 

Yridoet.-= 0.0: 

Nvldot = 0.0; 

% 

% Define constants of Nomoto’s first order equation 

% 

K = (Nvl*Yd1i - Yvl1l*Ndl1l)/[Yvl1* (Nrl-m*Xg) - Nvl* (Yrl-m) ]; 

Kv = -[((Nri-m*Xg)*Ydl - (Yrl-m)*Nd1]/{Yvl*(Nri-m*Xg) - Nvl*(Yrl-m)]; 
Tl = [(Yvldot-m)* (Nri-m*Xg) + Yvl* (Nrldot-Izz)]/[Yv1* (Nri-m*xXg) - Nvl* (Yri-m)}; 
T2 = [(Nvldot-m*Xg)*{(Yrl-m) - Nvl* (Yrldot-m*Xg)J]/[Yvl* (Nri-m*Xg) - Nvl*(Yrl-m)]; 
T3 = [Yd1* (Nvldot-m*Xg) - Nd1* (Yvldot-m)]/{Nvil*ydl - Yvl1*Nd1]; 
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T4 = [{(Nrldot - Izz)*Ydl - (Yridot - m*Xg)*Nd1)/{(Nrl-m*Xg)*Ydl1 - {(Yrl~m) *Ndl1]; 
T = T1+T2-T3; 
Tv = T1+T2-T4; 





% 

% Define storage vectors for Nomoto’s Model 

% 

kmaxl = (time/dt)+ 1; % Number of data points 

rl = zeros(1,kmaxl); % Initialize angular velocity vector 

rldot = zeros(1,kmaxl); %$ Initialize angular acceleration vector 

psil = zeros(1,kmaxl1); % Initialize heading angle vector 

psildot = zeros(1,kmaxl); %$ Initialize Yaw rate vector 

XL = zeros(1,kmaxl); $ Initialize X position vector (including side slip) 
XLdot = zeros(1,kmaxl1); % Initialize forward velocity vector(including side sli 
Pp) 

YL = zeros(1,kmaxl) ; % Initialize Y position vector (including side slip) 
YLdot = zeros(1,kmaxl1); %$ Initialize side slip velocity vector{including side s 
lip) . 
% Define storage vectors for Side Slip Model 

% 

v2 = zeros(1,kmaxl); % Initialize side slip veloctiy vector 

v2dot = zeros(1,kmaxl); % Initialize side slip acceleration vector 

X2 = zeros(1,kmaxl); % Initialize X position vector(including side slip) 
XZ2dot = zeros(1,kmaxl1); % Initialize forward velocity vector(including side sli 
Pp) 

Y2 = zeros(1,kmaxl) ; %* Initialize Y position vector (including side slip) 
YZdot = zeros(1,kmaxl); % Initialize side slip velocity vector{including side s 
Dap) 

% 

% Define storage vectors for Nomoto's Model w/correction factor 

% 


r3 = zeros(1,kmaxl) ; % Initialize angular velocity vector 

r3dot = zeros(1,kmax1l); $ Initialize angular acceleration vector 

psi3 = zeros(1,kmaxl); % Initialize heading angle vector 

psi3dot = zeros(1,kmaxl); % Initialize Yaw rate vector 

X3 = zeros(1,kmaxl) ; % Initialize X position vector(including side slip) 
% 


X3dot = zeros(1,kmaxl1); Initialize forward velocity vector(including side sli 
p) 


Y¥3 = zeros(1,kmaxl) ; % Initialize Y position vector (including side slip) 
Y3dot = zeros(1,kmax1); % Initialize side slip velocity vector(including side s 
lip) 

Timel = zeros({l,kmaxi); % Initialize time vector 

% 

for i = 1:kmaxl-l1; 

% 

$ Calculation loop for Nomoto’s Model 

% 

FLGOE (2). = (K*d = -r1(1)97T? 

r1(it+l)= rl(i) + rldot(i)*dt; 

psildot(i) = r1{i); 

psil(i+l)= psil(i) + psildot(i)*dt; 

XLdot(i) = cos(psil(i)); 


XL(i+1) = XL(i) + XLdot(i)*dt; 
YLdot(i) = sin(psil({i)); 
YL(i+1l) = YL(i) + YLdot(i)*dt; 


% 

% Calculation loop for Side Slip Model 

9 

vedot(i) = (Kv*d ~ v2(i))/TVv; 

v2(it+tl) = v2(i) + v2dot(i)*dt; 

X2dot(i) = cos(psil(i))- v2(i)*sin(psil(i)); 
X2(1+1) = X2(i) + X2dot(i)*dat; 

Y2dot(i) = sin(psil(i))+ v2(i)*cos(psil(i)); 
Y2({1+1) = Y2(i) + Y2dot(i)*dt; 

% 

* Calculation loop for Nomoto’s Model w/correction factor 
% 


R=0.5424*d+1.2560; 
r3dot(i) = (R*K*d - r3(i))/T; 
r3(i+1l)= r3(i) + r3dot(i)*dt; 


psi3dot(i) = r3(i); 

pSi3(i+l1)= psi3(i) + psi3dot(i)*dt; 
X3dot(i) = cos(psi3(i)); 

X3(1+1) = X3(i) + X3dot(i)*dt; 
Y3dot(i) = sin(psi3(i)); 


Y¥3(1+1) = Y3(i) + Y3dot(i)*dt; 
Timel(i+l) = Timel(i) + at; 
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% 


end 

% 

% NONLINEAR MODEL 

% 

%* Define hydrodynamic coefficients 

% 

Yous +0.5* PI*k + 2.4*s*Cb: 

YRS mi? = <b5"*s*Cch: 

Ybb = 1.25*(k/s) *(1-Cb) + 0.5; 

Yor = 0.277 (kss)*Cb = 0.07: 

VYorr “= .2.97"* (k/s)* (1 =cb):> 

VYbbr = 0.75" (ks) Cb. = 10.65: 

Nb = -k; 

NR = -0.54*k + k*%2; 

Nbb = -0.48%*(k/s)*(1-Cb) + .066; 

Nir 2:05" 6* Cb: = 09; 

Nbrr = -0.25*(k/s)*Cb + 0.05; 

Nbbr = =S7.5* (s*Cb) "2. # 138..4%s*Cb.— 1.6; 

Xudot = m-ml; 

Nrdot = Izz-I; 

Yvdot = m-m2; 

Yraot =" 0.03 

NVGGt:. =7505.°0% 

xo = 0.03 

Au. = 0% Os 

Oak Oe 

UR = Ue % effective rudder inflow speed 

% 

% Define storage vectors in advance to optimize performance 

% 

kmax = (time/dt)+1; %* Number of data points 

v = zeros(1,kmax); %* Initialize side slip veloctiy vector 
r = zeros(1,kmax) ; % Initialize angular velocity vector 
rdot = zeros(1,kmax); % Initialize angular acceleration vector 
b = zeros(i,kmax) ; % Initialize drift angle vector 

bdot = zeros(l1,kmax); % Initialize drift acceleration vector 
X = zeros(1,kmax); % Initialize X position vector 

Xdot = zeros(1,kmax); % Initialize forward velocity vector 
Y= zeros (1. kmax): % Initialize Y position vector 

Ydot = zeros(1,kmax) ; % Initialize side slip velocity vector 
psi = zeros(1,kmax); % Initialize heading angle vector 
psidot = zeros(1,kmax) ; % Initialize yaw rate vector 

Time = zeros({1,kmax) ; % Initialize time vector 

% 

%* Define nonlinear equations [calculation loop] 

$ 

for 1 = 1:kmax-1; 

Fn = ~(AR/(L*D))*Cn*UR*2*sin(d - gamma* (b(i)-(2*xR*r(i)))); %* Force acting on rudd: 
Xr = -(l-tr)*Fn*sin(q); 

Yr = -(1+aH) *Fn*cos (qa); 

Nr = -(xXR + aH*xH) *Fn*cos(d); 

% 

Xh = Xbr*r(i)*sin(b(i)) + Xuu *(cos(b(i)))%2; 

Yh = b(i)*Yb + r(i)*YR + xr(i)* (abs(r(i)))*Yrr + b(i)* (abs(b(i)))*Ybb + (b(i))*2*r(i)*Ybbr + b({ 


Pete CLO 2 Py bres 
Nh = b(i)*Nb + r(i)*NR + r(i)* (abs(r(i)))*Nrr + b(i)* (abs(b(i)))*Nbb + (b(i))*2*r(i)*Nbbr + b(3 
eam an Oo dead tl oh on oes 


% 

Yp = Q.0; 

Np = 0.0; 

Jp = (U*cos(b(1i))*(1l-wp))/(n*Dp); 

C1=0..523 

c2=-0.4861; 

e3=0 012127 

Kt = cl + c2*Jp + c3*dp%2; 

Xp = (Ctp* (1=-tp) *n*2*Dp*4*Kt) /(0.5*L*D*0%2).; 

% 

Xl = Xh + Xr + Xp; 

Yl = Yh + Yr + Yp; 

Nl = Nh + Nr + Np; 

% 

fl = (m-Yvdot)*(Izz-Nrdot) - (Nvdot-m*Xqg) * (Yrdot-m*Xg) ; 
Se Gi a «(ik oa se 

£3.= Nl= m*Xo*r(1); 

% 

rdot(i) = [(m-Yvdot)*f£3 + (Nvdot-m*Xg)*f£2)/f1; 
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r(itl) = r(i) + rdot(i)*dt; 

bdot(i) = [(Izz-Nrdot)*f£2 + (Yrdot-m*Xg)*f£3)/f1; 
b(i+tl)= b(i) + bdot(1)*dt; 

v(i) = U*tan(b(i)); 

psidot(i) = r(i); 


psi(itl) = psi(i) + psidot(i)*dt; 

Xdot(i) = cos{psi(i)}) - v(i)*sin(psi(i)); 
X(it1l}) = X(i) + Xdot({i)*dt; 

Ydot(i) = sin(psi(i)) + v(i)*cos(psi(i)); 
Y(itl) = Y(i) + Ydot(1)*dt; 

Time(i+1) = Time(i) + dt; 

end 

% 

* Plot results 

% 

figure(1) | 

DLlLOUCY X27 32" AL kby Hey yard? 
xlabel(’Transfer'),ylabel(’Advance’); 
gtext(’-- Nonlinear Model’); 

gtext(’. Nomoto model w/correction') ; 
gtext(’- Nomoto Model’); 

% 

figure (2) 
plot (Time (1:101) (12101) ;’=-" ,Timel (1:101) ,r1(1:101) , =", Timel (1:101) ,r3(22101).°.") grid; 
xlabel(‘Time(sec)‘),ylabel(‘Yaw Rate’) 
gtext(’-- Nonlinear Model’); 

gtext(’. Nomoto model w/correction’); 
gtext(’- Nomoto Model’); 

% 








4] 
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